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2 ; ABSTRACT 

, We discuss the ability of the planned Euclid mission to detect deviations from General Rela- 

■ tivity using its extensive redshift survey of more than 50 Million galaxies. Constraints on the 

gravity theory are placed measuring the growth rate of structure within 14 redshift bins be- 
tween z — 0.7 and z = 2. The growth rate is measured from redshift-space distortions, i.e. the 
anisotropy of the clustering pattern induced by coherent peculiar motions. This is performed 
in the overall context of the Euclid spectroscopic survey, which will simultaneously measure 
the expansion history of the universe, using the power spectrum and its baryonic features as 

' a standard ruler, accounting for the relative degeneracies of expansion and growth parame- 

\Q , ters. The resulting expected errors on the growth rate in the different redshift bins, expressed 

through the quantity /erg, range between 1.3% and 4.4%. We discuss the optimisation of the 
survey configuration and investigate the important dependence on the growth parameterisa- 
, tion and the assumed cosmological model. We show how a specific parameterisation could 

actually drive the design towards artificially restricted regions of the parameter space. Finally, 
in the framework of the popular "7 parameterisation", we show that the Euclid spectroscopic 
survey alone will already be able to provide substantial evidence (in Bayesian terms) if the 
growth index differs from the GR value 7 = 0.55 by at least ~ 0.13. This will combine 
with the comparable inference power provided by the Euclid weak lensing survey, resulting 
d i in Euclid's unique ability to provide a decisive test of modified gravity. 

Key words: cosmology - dark energy - large-scale structure of universe. 
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1 INTRODUCTION 

The current standard cosmological model, concordantly supported 
by virtually all available observations, tells us that we live in a low- 
density, expanding universe with a spatially flat geometry, that ap- 
pears to have recently entered a phase of accelerated expansion. 
The data require an extra mass-energy contribution in the form of 
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a fluid with equation of state w ~ — 1. This corresponds to adding 
in the equations of General Relativity (GR) a cosmological con- 
stant A, i.e. the term originally introduced by Einstein to obtain a 
static solution, thus building the standard Lambda Cold Dark Mat- 
ter model (LCDM). A cosmological constant, however, has a few 
disturbing features. The first feature is the fine tuning necessary 
to obtain its measured density value, which can be interpreted as 
the energy of the vacuum, that is extremely small compared to the 
corresponding scales of particle physics; the second issue is the So- 
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called coincidence problem: why, despite very different time evo- 
lutions, do dust-like matter and cosmological constant show com- 
parable densities today? Among possible solutions, scenarios with 
evolving 'dark energy' density from a cosm ological scalar field, 
have been proposed (see iFrieman et al, l l2008l for a review). Alter- 
natively, however, observations could simply indicate that it is the 
theory of gravity that needs to be revised (see lCopeland et al1l2006l 
for a review). 

These two radically different explanations cannot be distin- 
guished by measuring only the expansion history of the universe 
represented by the Hubble function H(z). A way to break this de- 
generacy is to look at the linear growth rate of density perturba- 
tions. This can be expressed as / = dlnG/dlna, where G(t) is 
the ti me-depende nt part of the solution of the linear growth equa- 
tion dHeatbll 19771) and a is the cosmic scale factor. Models with 
the same expansion history H(z), but base d on a different gravity 
theory predict a diffe rent growth rate f(z)_ (Maartens 2007; Linded 
l2005| - |Polarskill2006h (although see also lKunz & Sapondl2007T ~for 
possible issues). 

Measurements of galaxy clustering from large redshift sur- 
veys, quantified through the galaxy-galaxy correlation function (or 
its Fourier transform, the power spectrum P(k)), contain direct in- 
formation on both H(z) and f(z). 

Baryonic Acoustic Oscillations (BAO) within the last- 
scattering surface give rise to a characteristic feature in the 
galaxy distribution at comoving separations of ~ 150 Mpc (BAO 
have now been conc l usively seen in the clu stering of galaxies, 
e.g. ICole et alj 120051 ; lEisenstein etal] 120051 and more recently 
iKazin et alj |201C|) . This corresponds to the characteristic scale 
fixed by the comoving sound horizon at the drag epoch (shortly 
after recombination) and is accurately measured by Cosmic Mi- 
crowave Background (CMB) anisotropies (see recent measure- 
ments bv lKomatsu et alj|2009l ). Compared to the observed galaxy 
BAO peak position at different redshifts, this yields a measurement 
of H(z) in the radial direction and of the comoving angular diam- 
eter distance Da(z) in the transverse direction. 

Galaxy clustering as measured in redshift space also contains 
the imprint of the linear growth rate of structure, in the form of a 
measurable anisotropy. Such redshift- space distortion (RSD) is due 
to coherent flows of matter from low to high densities, linked to 
the growth of structure. When redshifts are used to measure galaxy 
distances, the contribution from peculiar velocities produces a dis- 
tortion of the clustering pattern, which at l inear scales is propor- 
tional to f(z) dKaiseJl987l : lHamiltonll 998). This can be measured 
by modelling the anisotropy of redshift-space two-point statistics, 
defining a proper correlation function £(r p , it) (or power spectrum 
P(fe|l , k±)), where r p and k\\ (it and k±) are the components paral- 
lel (perpendicular) to the line of sight. The anisotropy is caused by 
an additive term proportional to the variance of the velocity field, 
which can be parameterised by f(z)ag(z) (where erg is the rms 
amplitude of galaxy clustering) and provides an excellent discrim- 
inator of cosmological m odels, particularly if ag is normalised, for 
example using the CMB dSong & Percival 20091) . 

RSD were classically seen as a method to measure f?o (see ref - 
erence s in Hamilton ( 1998) or more r e cently Hawkins et al. (2003), 
iRoss et all d2007l), Ida Angela et alj d2008l) . ICabre & Gaztanagl 
< |2009|) . iDrinkwater et alj j2010l) ). In the dark energy context, they 
were initially seen as simply an effect to be correcte d for to ex- 
tract the full BA O information dSeo & Eisensteinl2003l> . Still in the 
context of GR, lAmendola et alj (2005) showed that the informa- 
tion contained in the GR growth function f (z) could improve er- 
rors on w(z) parameters by ~ 30% (see also Sa pone & AmendoTal 



120071) . As pointed out by iGuzzo et al] d2008l) , however, if no as- 
sumption is made on the gravity theory, RSD in fact provide us 
with a powerful test to test the dark energy vs modified grav- 
ity alternative by tracing the growth rate back in time. This was 
particularly interesting in the context of planned dar k energy sur- 
veys, and stimulated new interest in this technique dWang|l2008l: 
Lindei 2008; Nesseris & Pciivolaropoulos 2008; Acquaviva ct al. 
2008l ; ISong & Percivafe009l;IWhite et al.ll2008al ; |Percival & White 
2009). Following IGuzzo et alj 12008), RSD were suggested as a 
primary probe in the SPACE satellite proposal (the forerunner of 
what is now the Euclid spectroscopic probe, described in more de- 
tail be low) presented in response to the 2007 ESA Cosmic Vision 
Call dCimatti et alj|2009h . 

Currently on going s pectroscopic surve ys such as WiggleZ 
jBlake et alj20lTh , BOSS dWhite et alj20Tlh and VIPERS (Guzzo 
et al. 2012, in preparation), are mapping significant volumes of 
the distant universe and will produce measurements of fag cov- 
ering the redshift range out to z ~ 1.3 . Even larger and deeper 
galaxy surveys are planned for the next ten ye ars, both from ground 
and space. These are the proposed BigBOSS (Schlege l et al.|[2oTlh 
project using the refurbished KPNO and CTIO 4m telescopes, and, 
most importantly, the ultimat e redshift surve y from space by the ap- 
proved ESA mission Euclid dLaureiisl2 009) we already mentioned 
above. Main aim of these projects is precisely the solution of the 
dark energy puzzle. In particular, Euclid is a medium-size (M-class) 
mission of the ESA Cosmic Vision programme, which has recently 
been selected for implementation, with launch planned for 2019. 
Euclid will perform both a photometric survey in the visible and in 
three near-infrared bands, to measure weak gravitational lensing by 
imaging ~ 1.5 billion galaxies, plus a spectroscopic slitless survey 
of ~ 65, 000, 000 galaxies. Both surveys will be able to constrain 
both the expansion and growth histories of the universe. Their com- 
bination makes Euclid a unique experiment, with su perb precision 
and o ptimal cross-control of systematic effects (see lLaureiis et all 
1201 lh . In this paper we present the expected performances of the 
spectrosco pic survey, as desc ribed in the Euclid Definition Study 
ReporQ dLaureiis et al.ll201ll) , in quantifying the growth history 
and test for possible modifications of gravity as the origin of cosmic 
acceleration. 

We shall first explore the forecasted constraints on the growth 
rate both for the expected survey parameters, and for a pes- 
simistic galaxy density reduced by a factor of two. We will 
use the we ll-known parameterisation of the growth rate f(z) ~ 
\n m (z)r dPeebles|ll98d ; lFrvlll985l ; iLightman & Schechter|[l990l: 
Wang & Steinha rdJ 19981) , where Q. m = p m /p cr it, pm is the matter 
energy density and p cr it is the critical density making the universe 
spatially flat, to look at future constraints on the parameter 7, which 
characterises the gravity model, and test degeneracies arising when 
making assumptions on the background cosmology. 

We will then look at how constraints on growth are modified 
when survey specifications are changed, such as the total covered 
area and the corresponding galaxy number density (see discussion 
in Sec. [6} and the redshift range covered by the survey. In order to 
do this we need a measure of the quality of our growth constraint. 
For this reason we introduce a new simple Figure of Merit (FoM), 
analogous to t he dark energy FoM d efined in the Dark Energy Task 
Force Report dAlbrecht et al.ll2006h . We then show how this FoM 
depends on the already mentioned survey specifications. 

Expressing / as a function of 7 when considering sub-horizon 
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matter perturbation is a well justified choice for GR cosmologies 
and for some modified gravity theories (e.g. th e PGP model, a 
well-known five-dimensional model proposed in lDvali et al.ll2000l 
where the acceleration of the expansion is due to the leakage of 
gravity into an extra dimension) but it is not the only possibility. In 
some cases a co nstant 7 cannot describe the growth rate, as shown, 
for example, by |Pi Porto & Amendolal d2008h for a class of mod- 
els of dark energy coupled to matter. Therefore, when optimising a 
survey configuration for the growth rate measurement, it is impor- 
tant to consider not only one parameterisation of /, in order not to 
bias our decisions. The same was shown for the c ase of the equa- 
tion of state parameterisation in lWang et al.l d20ich . To have a more 
general perspective, we take a second growth parameterisation and 
compare the results obtained. We choose the param eterisation pro- 
posed bv lPogosian et "aH j2010h . lSong" et al.l d20 1 lh which is phys- 
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ical, works at any scale and it is also able to reproduce the DGP 
growth rate. 

Finally, we focus on the specific question "will the Euclid 
spectroscopic survey be able to distinguish between GR and mod- 
ified gravity models?" To answer it, we use a Bayesian approach, 
following the work of iHeavens e t al. (2007). Here, a way to fore- 
cast the Bayesian evidence is proposed, based on Fisher statistics. 
In particular, we look at how much 7 of a modified gravity model 
has to differ from the GR value in order for our fiducial survey to 
be able to distinguish it. 

Our work extends that o f lWangetal.l l l2010h , who considered 
using a Euclid-like survey to measure the dark energy equation 
of state, to furthe r consider the growth rate of st ructure. It con- 
firms t he results of ISimpson & Peacock! d2010h and lSamushia et ail 
d201 ll) . who specifically focused on cosmological model depen- 
dence. Parallel work t hat appeared while th is paper was in prepa- 
ration includes that of |Pi Porto et alj d201ll) . who look at bias and 
time-dependent parameterisations of 7, while we specialise on dif- 
ferent growth parameterisat i ons an d on the dependence on survey 
spe cificationslBelloso et alj d201 lh . who find an exact solution of 7 
and lZiaeepouj J2OI lh . who estimate constraints from wide surveys 
based on a new growth parameterisation for modified gravity and 
interacting dark energy models. A difference between our work and 
all the papers mentioned ab ove is that our resul ts are updated to the 
latest Euclid configuration dLaureiis et al.l201 lb , which includes re- 
sults from the most accurate simulations of the instrument. 

The outline of the paper is the following. In Sec. [2] we de- 
scribe RSP as an observable to measure growth, while in Sec.[3]fhe 
Fisher matrix method used to forecast errors on growth and other 
cosmological parameters is outlined, and the fiducial cosmology 
is defined. Sec. [4] describes how the survey was modelled, while 
Sec. [5] shows our forecasts, and in the following Sec. [6] we define 
the growth Figure of Merit, which we apply to our forecasted sur- 
vey data. In Sec. [7] we compute the forecasted Bayesian evidence 
and we finally conclude in Sec. [8] Appendix[A]contains a brief re- 
view of the theory of the growth of linear small-scale perturbations 
and of growth parameterisations used in the literature. 



2 MEASURING GROWTH: REDSHIFT SPACE 
DISTORTIONS 

Our observable, containing both BAO and RSP, is the galaxy 
power spectrum. We write the observed power spectrum P b s as 
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where the subscript r indicates quantities in the reference 
cosmology chosen to compute the power spectrum, = 
k- L D A {z)/D A (z) I , fcj = k i H(z) I /H(z) are the wave 
modes perpendicular and parallel to the line-of-sight, k — 
V k ±2 + fcli 2 , /1 = k ■ r/k = k^/k and P s hot is a scale- 
independent offset due to imperfect removal of shot-noise. The 
term [D A (z) T /D A (z)] 2 [H(z)/H(z) T ] represents the distortion 
of the power spectrum due to the Alcock-Paczynski effect 
dAlcock & Paczvnski|[l979]) . which also has the impact of chang- 
ing the true (k, fi) into reference (k r , fj, T ). 

Notice here that the term [P ma ttcr(fc)/o"g] z _ is independent 
of the normalisation of the matter power spectrum. 

We allow for the growth of structure to vary away from that of 
a LCPM model in two directions. First we allow the growth rate, 
which is well appr oximated by / = QZi (/ = + (7 — 4/7)Qk 
for curved space, iGong et al.l l2009h with 7 = 0.545 in LCPM 
models, to have 7 7^ 0.545. Second, we directly allow a modifi- 
cation of the standard equation for the ti me-like metric potentia l 
ty. This modificat ion is described following Pogosian et alj d2010l) : 
ISong et alj d201 l|) by a constant parameter /i s 7^ 0. This choice 
and its motivation are discussed in more detail in Appendix A. We 
will use these parameters to model growth from observations of 
redshift-space distortions. 



3 FORECASTING THE ERRORS 

In this section we give details of the method used to estimate the 
expected minimum errors obtainable on 7 or fi s and the other cos- 
mological parameters from future redshift surveys, once the mea- 
surement error on the observables is known. 

To perform our forecasts we use the standard Fisher ma- 
trix approach which was introduced to f orecast errors on P(k) 
and derived parameters dTegmarkl [l997h and then adapted by 
ISeo & Eisensteinl d2003h to the measurement of distances us- 
ing the BAO position, avoiding the "noise" of RSP. A se- 
ries of papers in the literature have then used and discussed 

this technique to extract the information on both expansion and 

I 1 1 1 I — 1 1 — 1 

growth (for example rec ently Wang 2008; White et al. 20081 

Simpson & Peacock 20ld : IWang et al . 2010; Sa mushia et a l. 201 
|Pi Porto et alj|201 ll). Here we follow in particular the approxima 

tions discussed in Samushia et all d201lh . 

Briefly, the Fisher matrix is defined as (see e.g. iBassett et alj 

2009) Fij = - (d 2 In C/dpidpi) , where C is the likelihood func- 
tion and pi are the model parameters whose error one wishes to 
forecast. It is possible to show that the Fisher matrix is the inverse 
of the covariance matrix, Fij = , in the (strong) assumption 
that the likelihood is a Gaussian function of the parameters and not 
only of the data. The galaxy pow er spectrum Fisher matrix can be 
approximated as (Tegmark 1997) 
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where P b s (k) is the galaxy power spectrum, Eq. ([T}, and its 
derivatives are computed in a chosen fiducial model, V e ff = 
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Vsurv [1 + l/(nP bs(k))]~ 2 , Vsurv is the volume of the survey 
and n is the gal axy number density. 

Following lSamushia et al.l J201 ll) . we neglect P s hot. since this 
term should only introduce negligible error. We also neglect the 
dependence of P m atter(fc) on flbh 2 , £l m h 2 , h and n a . This means 
that we do not use the information coming from the shape of the 
matter power spectrum. We finally multiply the integrand of the 
Fisher Matrix by a factor accounting for a possible error in red- 
shift, exp[— (fc re f/x rc f ca z /H Te {) 2 )], where a z = 0.001(1 + z) is 
the standard deviation of the redshift error expected in Euclid. In 
each redshift bin centre d at z — Zi we compute a Fisher matrix 
dSeo & EisensteirJl2003l) with parameters pj 

Pj = {f( z i) a a( z i) > b(zi)a$(zi) , \nD A (zi) , InH(zi)} , (3) 

where pj in different redshift bins are considered to be independent; 
then all the matrices corresponding to all &i are summed. From the 
total Fisher matrix it is possible to estimate the errors on each pj 
in each redshift bin, by marginalising over all other parameters. 
To obtain errors on cosmological parameters, we marginalise over 
b(zi)as(xt) and project the obtained matrix into the final cosmo- 
logical parameter set. 

Regarding the latter, we test three nested models: 1) a simple 
quasi-LCDM with a constant w fixed to —0.95 instead of —1 as 
in LCDM, which we dub qLCD]V0, 2) a model where dark energy 
has constant equation of state w, dubbed wCDM and 3) a model 
where the equation of state of dark energy is allowed to vary, fol- 
lowing the evolution w = wp + w a (l — a), which we call CPL 
dChevallier & Polarskil200ll ; lLinderl2063h . For all these models we 
consider both the flat and the curved space cases. The full set of pa- 
rameters, from which we pick the appropriate subset according to 
the chosen model, is therefore 

q 3 = {h, fi m , fifc , wq , w a , 7 or Ms , °~a(z = 0)} , (4) 

where h = Hp/ (10 kms~ 1 Mpc~ 1 ). This is different from 
ISamushia et al.l d201lh since they do not consider erg (2 = 0) and 
their resulting errors will be consequently smaller with respect to 
this work, as we will see in Sec. [5] 

Our fiducial model is, as in the Euclid Definition Report, a 
flat consta nt w cosmology wit h parameter values as the best-fit 
WMAP-7 jKomatsu et alj201ll) results (except for w): 

h = 0.703 Q m = 0.271 ft fe = 0. a s (z = 0) = 0.809 

w Q = -0.95 w a = 0. 7 = 0.545 //x s = 

ft b = 0.045 n s = 0.966 (5) 

The fiducial ma tter power spectrum for this cosmology is computed 
using CAMB0 dLewis et al j2000h . 

We assume a scale-independent bias, which is a good approx- 
imation for large enoug h scales. As a fid ucial bias, we take the 
bias function derived by lOrsi et all J2009t> using a semi-analytical 
model of galaxy formation. For an analysis of how the bias can 
be constrained and of the impact of assuming a biasing model on 
the estimates of the growth factor we refer to the parallel work of 
|Pi Porto et a 1). Number densities and biasing parameters are 
summarised in Table [T] together with the integration limits in k. 
The latter correspond to scales R such that cr 2 (R) — 0.25, with an 
additional cut at fc max = 0.20 h Mpc~ x . 

2 We take w ^ —1 to match the choice of Laureiis et al. (2011), moti- 
vated by the possibility of computing cosmological perturbations avoiding 
the w = — 1 barrier. 

3 http://camb.info 



In a separate paper, Bianc hTet al.l | |2012|) perform a detailed 
study of statistical and systematic errors in RSD measurements, 
using a large set of mock surveys built from numerical simulations. 
They compare their results with predictions from the Fisher matrix 
code used here, finding fairly good agreement when considering 
only linear scales. 



4 MODELLING THE EUCLID SURVEY 

Euclid wi ll cover an area of 15 , 000 deg 2 in both imaging and spec- 
troscopy (Laureiis et al. 2011), measuring redshifts in the infrared 
band (0.9 — 2 /im) for ~ 65 Million galaxies, using a slitless spec- 
trograph. With this technique, the redshift measurement relies on 
the detection of emission lines in the galaxy spectra, which in the 
chosen wavelength range in the vast majority of cases will be the 
Ha line, redshifted to 0.7 < z < 2. 

Intrinsic to the slitless technique is the impossibility to de- 
fine a priori a survey flux limit (as normally done in a classical 
slit survey, which is based on a well defined target sample selected 
to a given magnitude or flux limit). Spectra are intrinsically con- 
fused by superpositions and the actual flux limit depends not only 
on the nominal signal-to-noise reachable through a given exposure, 
but also and fundamentally on the strategy devised as to resolve 
the confusion among the different spectra. In the case of Euclid, 
this is achieved first by splitting the wavelength range into two sub- 
exposures through a "blue" and a "red" grisms, covering respec- 
tively the wavelength ranges 0.9 — 1.4 /im and 1.4 — 2.0 fim; this 
has the advantage of halving the length of the spectra on the de- 
tector, thus reducing superpositions. Secondly, the blue and red ex- 
posures are in turn split into two sub-exposures, observed rotating 
the field of view by 90°. This makes for four different exposures of 
a given field, which result in a sensible treatment of spectral con- 
fusion. To verify this and compute the success rate of the survey 
(i.e. the fraction of correctly measured redshifts over the total num- 
ber of spectra), it was necessary to develop since the early stages 
of the project an advanced end-to-end simulation pipeline. An ac- 
curate description of these simulations is beyond the scope of this 
paper, and we refer th e reader to the Euclid Definition Study Report 
dLaureiis et al.ll20lTh . and to the specific paper (Garilli et al. 2012, 
in preparation). Here we recall only the main concepts, which are 
relevant for the present analysis. 

The end-to-end spectroscopic simulations take into account 
Euclid's instrumental (point-spread function, resolution and instru- 
mental background) and observational (exposure time, astrophysi- 
cal background) parameters, to build an artificial "observation" of 
a realistic galaxy field. The input data set is buil t starting from th e 
COSMOS catalogue of photometric redshifts dllbert et al] 2009), 
which contains all relevant information (coordinates, redshift, lu- 
minosity, spectral energy distribution, stellar formation rate, etc.). 
Most importantly, the COSMOS depth and extensive wavelength 
coverage (over 30 spectral bands), allows us to assign a well- 
defined spectral type and thus a realistic distribution of Ha equiva- 
lent widths, down to large distances. It is also important that, being 
based on real observations, the catalogue includes also a realistic 
clustering of the sources. Stars are then added upon the galaxy cat- 
alogue. The package aXeSIM3 is then used to generate the four 2D 
dispersed images of the field and then extract the four ID spectra 
for each target. Redshift are measured making use of RESS (Red- 
shift Evaluation from Slitless Spectroscopy), an automatic software 
running within IRAF environment, which has been devised (Ros- 
setti E. 2012, in preparation) to reduce and analyse highly contam- 
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Figure 1. Predicted mean number density of galaxies in each redshift bin 
centred in z, expected from the baseline Euclid wide spectroscopic sur- 
vey, given the instrumental and survey configurations and the estimated ef- 
ficiency. 



inated slitless spectra. RESS' z measure is currently implemented 
only for low z galaxies (0.7 < z < 2.0) and its extension for high z 
objects is in progress. Redshift evaluation is based on the position 
of the Ha line and any other emission lines, when detected, for 
which a flux is also measured. A reliability flag for each measured 
redshift is then obtained by further pr ocessing the spec tra through 
the EZ redshift measurement code dGarilli et alJbOld) . Compari- 
son of the input and output catalogues allows one to estimate the 
success rate of the survey in terms of completeness and purity as 
a function of redshift and Ha flux (see Euclid Definition Study 
Report, Fig. 6.10). Rather than trusting the absolute redshift dis- 
tribution emerging from the simulated field, a more conservative 
choice is to use this output as weight, to be applied to the most 
up-to-date predict ions for the redshift distribution of Ha emitters 
dGeach et alj|2008l) . This produces the expected distribution of the 
number of galaxies with measured redshift in each redshift bin. 
From this one can calculate the galaxy number density at each z, 
which is shown in Fig.Q]for our fiducial cosmology of Eq. ((5}. 
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Table 1. Galaxy biasing parameter b and fc max of integration for each red- 
shift bin centred in z for the Euclid spectroscopic survey baseline configu- 
ration, having an observed area of 15, 000 deg 2 



The values of / ag are com puted in the case of lGuzzo et al. I d2008l) 
and lHawkins et al] d2003l) by using the value of f /b given by the 
authors and computing bag from b and th e reference cosm ology 
they ad opt for the computati on of b (or of lLahav et al.|2002l i n the 
case of lHawkins et alJl2003h : in the case of iRoss et alj ( 120070 bag 
was computed using the expressiorfj] ( Zehavi et alj200^). (bag) 2 = 
So d V V 2 (3 - 9y/4 + 3v 3 /16). ICabre & Gaztanagal Jiooj 

indicate directly their value of bag, while IB lake et al] ( 1201 lh and 
ISamushia et al] d201 ll) compute directly fas. Error bars are ob- 
tained through the error propagation formula for uncorrected data, 
when not directly specified in the papers. 

Together with the (solid black) curve representing our fidu- 
cial f ag, we also show for comparison a (dashed green) line for 
flat DGP, (calculated by numerical integration of the correspond- 
in g equation for /) and a (dotted red) line for the coupled model 
of Di Porto et al. (2011), computed using the parameterisation of 
|Pi Porto & Amendolal d2008l) with a coupling f3 c = 0.2 (both with 

= 0.271 and the same ag(zcMB) of our fiducial model). 

We notice that we reach accuracies between 1.3% and 4.4% 
in the measurement of / ag depending on the redshift bin, where 
the highest precision is reached for redshifts z ~ 1.0. 



5 STANDARD PREDICTIONS FOR EUCLID 

For our computations here, we split the Euclid predicted redshift 
distribution over the range 0.7 < z < 2, into 14 bins with 
Az — 0.1. Using the predicted galaxy number density in each bin 
shown in Fig.Q] we obtain the error on our observable, the power 
spectrum, and estimate the resulting precision on the measurement 
of / ag after marginalisation over the other parameters. We plot er- 
rors on / ag in Fig. [2] (dark blue error bars), where we also show 
for comparison current measurements of / ag (light pink and ma- 
genta error bars) and the pessimist ic case of observin g only half 
the number of galaxies forecasted in lGeach et al.ld2008l) (light blue 
error bars), as the authors themselves claim that their counts may 
be wrong by a factor of 2. 

Current measurements shown in Fig. [2] are listed in Table [2] 



5.1 Comparison to other surveys 

Together with Euclid, other ongoing and future surveys will con- 
strain cosmology by measuring fag. Here we compare the rela- 
tive errors on fag obtained using different spectroscopic galaxy 
redshift surveys. In p articular, we consider the BOSS surve)[j (see 
ISchlegel et al]|2009l) and the BigBOSS3 Emission Line Galaxies 



4 This formula actually gives us the non-linear bag, since we have used 
the non linear estimate of £ of Ross et al. (2007) to compute it. What we 
needed to obtain the linear fag would be the linear bag, bu t we do not 
have it. Therefore our estimate of fag for the lRoss et alj 120071) datapoint 
might be 5 — 10% higher than it should. 

5 http://cosmology.lbl.gov/BOSS/ 
s http://bigboss.lbl.gov/ 
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Figure 2. Fisher matrix forecasts of the errors expected on the growth rate (dark-blue error bars), expressed through the bias-free combination f(z)ag(z), 
obtainable from the Euclid redshift survey through the combination of amplitude and redshift-space anisotropy of galaxy clustering. The light-blue error bars 
(shown with a slight offset in redshift for visualisation purposes) represent the case of a galaxy density reduced by a factor of two with respect to that forecasted 
for the galaxies observed by Euclid (Geach et al. 2008). The solid black line represents the fiducial / as, computed for the cosmology shown in Eq. {5}. The 
dashed green line shows the growth of a flat DGP model (calculated by numerical integration of the corresponding equation for f(z)). The red dotted line 
represents / as of a coupled model with coupling parameter C = 0.2. All models are computed for Q m = 0.271 and for the same ag(zcM b) as for the 
fiducial model. In the same plot we also show measurements of / as from past surveys (magenta error bars) and the recent WiggleZ survey (pink error bars), 
see explanation in the text. 



survey 


reference paper 


z 


f&8 


WDS F22 
wide 


Guzzo et al. (2008) 


0.77 


0.49 ±0.19 


2SLAQ 
galaxy 


Ross et al. (2007) 


0.55 


0.50 ±0.07 


SDSS LRG 


Cabre & Gaztanasa (2009) 
Samushia et al. (2011) 


0.34 
0.25 
0.37 


0.53 ±0.07 
0.35 ±0.06 
0.46 ± 0.04 


2dFGRS 


Hawkins et al. (2003) 


0.15 


0.39 ±0.08 


WiggleZ 


Blake et al. (2011) 


0.22 
0.41 
0.6 
0.78 


0.49 ±0.07 
0.45 ± 0.04 
0.43 ± 0.04 
0.78 ± 0.04 



Table 2. Current measurements of fas 



(ELGs) and Luminous Red Galaxies (LRGsQ- Regarding the fidu- 



7 We thank the BigBOSS consortium for providing their latest yet unpub- 
lished estimate of their expected galaxy densities, which we used in creating 
this plot. 
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Figure 3. Relative error on / erg of Euclid (dark-green circles, light-green 
circles for the pessimistic case of half the galaxy number density), BOSS 
(dark-red squares), BigBOSS ELGs (blue triangles) and LRGs (orange dia- 
monds). 



cial bias, we use the forecasts by lOrsi et al.l d2009h for BigBOSS 
ELGs. We use b = 2G(0)/G(z) (where G{z) is th e standard 
linear growth rate) for BOSS and BigBOSS LRGs (see iReid et al.l 
feOlCh ). Table[3] summarises the main characteristics of these sur- 
veys. 
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survey 


redshift range 


area [deg 2 ] 


n [h 3 Mpc" 3 ] 


bias 


BOSS LRG 


0.05 < 2 < 0.65 


10, 000 


3 x 10~ 4 


2.0G(0)/G(z) 


BigBOSS LRG 


O.K z < 1.1 


14, 000 


unpublished (see footnote 7) 


2.0G(0)/G(z) 


BigBOSS ELG 


O.K z < 1.8 


14, 000 


unpublished (see footnote 7) 


see Orsi et al. (2009) 



Table 3. Future and ongoing galaxy redshift surveys and their main properties 
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Figure 4. Forecasted errors on cosmological parameters using the Euclid spectroscopic survey. Solid (black) lines correspond to qLCDM, dashed (red) lines 
to wCDM and dot-dashed (cyan) lines to CPL. Thick lines represent models with = 0, while thin lines represent curved models. Left plot: marginalised 
errors on Q m and 7. Right plot: marginalised errors on Wq and 7. 



The results are shown in Fig. [5] We first notice that Euclid 
(represented by dark-green circles) will obtain the most precise 
measurements of growth, even in the pessimistic situation of detect- 
ing only half the galaxies (light-green circles). In redshift coverage 
it will be perfectly complementary to BOSS. The partial overlap 
with BigBOSS, whose ELG sample will reach similar errors up to 
z ~ 1.4, will allow for interesting and useful independent mea- 
surements and cross-checks. 



5.2 Cosmological parameters 

We next look at how the errors on / erg, H and Da project into 
errors on cosmological parameters. We only show results for 7 
and briefly mention results for [i s (see Appendix [B] for more de- 
tailed error forecasts), while in the next section we will compare the 
two parameterisations 7 and /j, s . In order to understand the model- 
dependence of these forecasts, we use the nested cosmologies de- 
scribed in Sec. [2] the more complicated model is a generalisation 
of the less complicated one, and the generalisation consists in the 
addition of one extra parameter. We also study the influence in the 
estimate of background parameters of assuming GR or allowing for 
a different constant 7. Results are shown in Figs. 14171 

From Fig. [4] and [5] and by comparing left with right panel of 
Fig.|6]it is clear that what affects most strongly the error estimate, 



by enlarging the ellipses and even changing the degeneracy direc- 
tion of parameters, is the assumption of different dark energy mod- 
els. The assumption of zero curvature also has an influence. The 
less parameters the cosmological model possesses, the stronger is 
the influence of fixing Qk, as can be seen from Fig. [4] (compare 
thick and thin lines having the same line style/colour) and Fig. [6] 
(compare solid and dashed lines with same thickness). = 
affects more strongly the measure of background parameters, and 
only indirectly constraints on 7 (note the thickening of error el- 
lipses in Fig.|4]in the direction of background parameters Q m and 
wo, respectively). 

As regards the /j, parameterisation, the error on fi s is fore- 
casted to be between ~ 0.5 and ~ 2, depending on the dark energy 
model. The assumptions on w and those on Qk affect more strongly 
constraints on this growth parameter, in the case of the /i s — fl m 
error ellipses. 

The assumption of GR has a different influence on different 
pairs of parameters. As regards Q m — fl^, we can see from Fig. [5] 
that fixing 7 reduces only the error for qLCDM, while errors on 
wCDM and CPL are unaffected. Errors on wo — w a are very weakly 
affected by the assumption of GR and also of zero curvature (see 
Fig. |7j, so that their determinati on is rather robust. Th is conclu- 
sion agrees with what obtained in Samushia et al. (see their 
Figs. 4 (a) and 5 (a) of the final published version), where the small 
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Figure 5. Forecasted marginalised errors on Q m and f2^ using the Eu- 
clid spectroscopic survey. Solid (black) lines correspond to qLCDM, dashed 
(red) lines to wCDM and dot-dashed (cyan) lines to CPL. Thick lines corre- 
spond to models where 7 = 0.545 (its GR value) while thin lines represent 
models where 7 is allowed to assume other (constant) values. Note that for 
wCDM and CPL the thick and thin lines are one on top of each other: fixing 
7 does not have any impact on constraints on f2 m and f!fc. 



differences are due to the presence of an extra parameter in our 
analysis over which we marginalise, namely as(z = 0) and to the 
different survey specifications. The different degeneration direction 
of Q, m — Qk is due to the wrong sign convention being adopted in 
figure 5 of Samus hia et al. 1) which had the effect of inverting 
the ilfc a xis. We have checked t hat our code gives exactly the same 
results as Samus hia et al.l d201ll) if the parameter erg (z = 0) is fixed 
and the survey specifications are identical. The pair of parameters 
which is most affected by the choice of fixing 7 is fl m — Wo, more 
strongly in the case of wCDM, as can be seen in Fig. [6] 

5.3 Adding Planck 

We then quantify the impact of adding Planck constraints to the 
Euclid spectroscopic survey constraints, to check how much errors 
reduce. The Planck satellite (part of the Cosmic Vision programme 
by the European Spac e Agency) has bee n launched in 2009 and it 
is presently operating (Taub er et al.l201oT) . Its results on cosmology 
will be made public in the next years and will be the state-of-the-art 
data on CMB temperature and polarisation. It is natural and conve- 
nient to combine these data to Euclid galaxy survey data, first of 
all because these probes are highly complementary, observing the 
early and late universe respectively, and secondly because of the 
presen ce of the BAO feature in both datasets. As in lSamushia et al.l 
(2011), we utilise the Dark Energy Task Force Planck Fisher, com- 
puted for 8 parameters: h, f2 m , fij., wo, w a , (Ts{z — 0), n B and 
flb. Since this Fisher matrix is computed assuming GR, to gener- 
alise it for arbitrary 7 in order to sum i t to our galaxy survey Fisher 
matrix, we use the same method as in Samushia et al.l d201ll) (see 
their Appendix B for details). In Figs. 1811 Fl we compare errors with 
(solid purple lines) and without (red-dashed or cyan dot-dashed 
lines) adding Planck, for a CPL dark energy. 

As can be seen from Fig. [8] Planck noticeably improves 
marginalised errors on 7, obviously by constraining the background 
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Figure 7. Forecasted errors on the CPL parameters wq and w a using the 
Euclid spectroscopic survey. Solid (black) lines correspond to models with 
Qk = 0, and dashed (red) lines to curved models. Thick lines refer to 
models where 7 = 0.545, while thin lines represent models where 7 is 
allowed to take other (constant) values. 



parameters. The error most reduced by Planck is, as one expects, 
that on Q m — (see Fig. [9), quite independently of the assump- 
tions on growth (compare thick lines-corresponding to assuming 
GR, to thin lines-where 7 is not fixed). As regards constraints on 
the alternative /1 parameterisation, here the impact of adding Planck 
data on combined /j, s — Q m and fi s — wo constraints is stronger than 
for 7: the projected error on /j, s is reduced from ~ 2 to ~ 0.7. This 
is shown in more detail in Appendix [B] loint Euclid-Planck con- 
straints on Q. m — wq (Fig. I lOt and wo — w a (Fig. II It depend more 
strongly on the assumption on 7, but are in any case a decisive 
improvement with respect to Euc lid-only constraints. Th is is again 
consistent with what obtained bv lSamushia et al.l d201 ll) . who find 
slightly tighter constraints for the reason explained above. 



6 A FIGURE OF MERIT FOR GROWTH 

Having estimated errors on parameters in the fiducial Euclid con- 
figuration, in this section we move to investigating how well the 
growth of structure is measured when varying essential survey pa- 
rameters, such as the observed area and the corresponding galaxy 
number densities observed, as well as the redshift range where 
galaxies are observed. To quantify the success in measuring growth, 
we resort to a dedicated Figure of Merit, analogous to the Dark 
Energy Task Force F oM. In the Dark Energy Task Force report 
dAlbrecht et alT 2006). a FoM was introduced to quantify progress 
in measuring the properties of dark energy. This was defined as the 
inverse of the area enclosed in the 95% confidence level contour of 
the equation of state parameters wo — W a - 

6.1 Definition of the FoM 

Here we introduce an analogous FoM aimed at quantifying the pre- 
cision of the growth measurement. Our FoM is defined as the in- 
verse of the area enclosed in the 95% c.l. contour of the growth 
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Figure 6. Forecasted errors on Q m and wo with the Euclid spectroscopic survey. Solid (black) lines correspond to models with = and dashed (red) 
lines to curved models. Thick lines refer to models where 7 is fixed to the GR value, while thin lines represent models where 7 is allowed to vary. Left panel: 
wCDM model. Right panel: CPL model 



0.70 - 




0.40 - 



0.20 0.25 0.30 0.35 



Figure 8. Comparison of marginalised errors on curved CPL, obtained using 
correspond to error ellipses using Euclid data only, while solid (purple) contours 
on Q m and 7. Right plot: marginalised errors on uiq and 7. 
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Euclid only and joint Euclid and Planck data. Dot-dashed (cyan) contours 
show joint constraints from Euclid and Planck. Left plot: marginalised errors 



parameter (7 or /i s ) and Q. m . We choose Q. m as our second param- 
eter because we want the FoM to be usable for all models, while 
e.g. the equation of state, which could represent an alternative to 
Q. m , is fixed to —0.95 in qLCDM and therefore cannot be used. 
Moreover, our main observable, the RSD, constrains Q m very ac- 
curately. The FoM as we define it can be computed through this 
simple formula: 



FoM = 



y/ det F-,n ri 



-2tt In 



(6) 



6.2 Dependence of the FoM on survey specifications 

We compare here 7 and ^t s and explore three survey parameters: the 
survey area, corresponding to an associated galaxy number density, 
the maximum redshift z max and the redshifts range z m in — z max - 

If we wanted to vary the survey area independently of the other 
survey parameters, we would find a simple proportionality relation: 
FoM oc area, independently of the parameterisation assumed. This 
can be easily understood by looking at Eq. ((2} and knowing that 
Vsurv is always proportional to the total survey area, which means 
that Fij oc area. Therefore, det F-y n m oc area 2 and from Eq. ©, 
FoM oc area. 



or through an analogous one for fi s . 



A more realistic question to be asked when planning a survey 
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Figure 9. Comparison of marginalised errors on background parameters 
f2 m and of a curved CPL model, obtained using Euclid and joint Euclid 
and Planck data. Dot-dashed (cyan) contours correspond to error ellipses 
using Euclid data only, while solid (purple) contours show joint constraints 
from Euclid and Planck. Thick lines refer to models where 7 = 0.545, 
while thin lines correspond to models where 7 is allowed to vary. Thin and 
thick cyan lines are superimposed because the CPL constraints from Euclid 
only data do not change when fixing 7. 
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Figure 10. Comparison of marginalised errors on background parameters 
Q m and u>o of a curved CPL model, obtained using Euclid and joint Euclid 
and Planck data. Dashed (red) contours correspond to error ellipses using 
Euclid data only, while solid (purple) contours show joint constraints from 
Euclid and Planck. Thick lines refer to models where 7 = 0.545, while 
thin lines correspond to models where 7 is allowed to vary. 



is rather whether it is more convenient to invest the total survey 
duration by mapping a smaller area for a longer time, so to increase 
the number density of observed galaxies, or rather mapping a larger 
region of the sky but collecting a sparser sample of galaxy spectra. 

To understand this we have built a simple linear scaling on 
how the number of galaxies per redshift bin increases when the 



Figure 11. Comparison of marginalised en'ors on the curved CPL back- 
ground parameters wo and w a , obtained using Euclid and joint Euclid and 
Planck data. Colours and line styles are the same as in Fig. |10| 



area is reduced and the exposure time per sky patch is increased but 
the total available time is fixed. In particular, the same end-to-end 
simulations described in Sec. [4] provide us also with an estimate of 
the number density increase obtained when increasing the observ- 
ing time per observed patch by 25%. We assume that this increase 
corresponds to a 25% decrease in survey area covered, given that 
the total survey time available is constant, and that this relation is 
linear. From these considerations, we derive the following formula: 



n(z) — 7115000(2) 



C(z)) + C{z) 



15,000 deg 2 



(7) 



where ni 5 ooo is the number density corresponding to an area of 
15, 000 square degrees and the values of C(z) are listed in Table 
[4] More precise answe rs would be obtained using e.g. the method 
of Bassett et al. ( 2005) but we expect this approximation to be suf- 
ficient at this stage to roughly describe the trade-off between area 
and number density. In Figure Q/2] we show the dependence of the 
FoMs for 7 and /i s (left and right panels respectively) on the area - 
galaxy number density. In order to better appreciate the functional 
behaviour, we plot the relative improvement in the FoM when 
varying the area with respect to the FoM of area = 12, OOOdeg 2 : 
(FoM(area) — FoM(area = 12000deg 2 )) /FoM( area = 
12000deg 2 ). Here we can see that the improvement in the growth 
measurement is nearly linear, at least for the area interval con- 
sidered (12, OOOdeg 2 - 20, OOOdeg 2 ), and mildly dependent on 
the dark energy model considered, but also rather slow, so that a 
large increase in area might not be favoured in case the cost for 
obtaining it grows too fast. Comparing left to right plot, we also 
note that in the case of the fi s parameterisation the improvement 
of the FoM depends more strongly on the cosmological model one 
has previously assumed. The largest improvement appears for flat 
CPL. 

As regards the dependence of the FoM on z max , we evalu- 
ated it by calculating our total Fisher matrix as the sum of the 
Fisher matrices computed for the redshift bins from z = 0.7 
to z — Zmax only (see Sec. [3] for further details on the gen- 
eral Fisher matrix calculation technique). Fig. 1131 shows the rel- 
ative difference of the VdM{z max ) to the FoM for z ma x = 1: 
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Figure 12. Dependence of the relative FoM on the survey area and corresponding galaxy number density. The FoM is computed for the growth parameters 7 
(left panel) and )j, s (right plot) and normalised to the FoM for area = 12000deg 2 : (FoM - FoM(area = 12000deg 2 )) /FoM(area = 12000deg 2 ). Blue 
circles, red squares and green triangles are for qLCDM, wCDM and CPL respectively. Full (empty) symbols represent flat (curved) models. Solid (dashed) lines 
join flat (curved) models. The arrow indicates the baseline area of the Euclid survey, 15, 000 deg 2 , resulting from an optimisation of the joint spectroscopic 
and photometric surveys. 



Table 4. Trade-off between number density and area: C(z) for each redshift 
bin, to be inserted in Eq. j7), in order to relate the increase in number density 
to variations in the total survey area. 



z 


C(z) 


0.7 


0.365786 


0.8 


0.400569 


0.9 


0.408442 


1.0 


0.412037 


1.1 


0.418330 


1.2 


0.418001 


1.3 


0.422914 


1.4 


0.422946 


1.5 


0.439024 


1.6 


0.490589 


1.7 


0.506262 


1.8 


0.521747 


1.9 


0.531328 


2.0 


0.551155 



(FoM(2: ma;c ) - FoM(z max = 1)) /FoM(z ma x = 1). Here the 
amplitude of the spread due to model assumptions is much wider 



than for the area - galaxy number density dependence, and wider in 
the case of 7. The latter is because the derivative of our observable 
/08 with respect to fjb s decreases with redshift while its derivative 
with respect to 7 increases (as we have checked numerically), so 
that at z ~ 1 fas is more sensitive to changes in 7 than in fi s . This 
explains why in the case of 7 and for the CPL model, which has 
the strongest time-dependence of the background cosmology, there 
is more advantage from higher redshift bins than for the case of fj, a - 
In both parameterisations complex models gain more from higher 
redshift data, so e.g. by increasing z„ lax wCDM's FoM improves 
more than qLCDM's does, and the FoM of models with 7^ 
improves more than the corresponding one for flat models (with the 
exception of CPL with 7). In the case of 7, this is easy to explain. 
We know that / = f2 m (z) 7 . In the case of a flat qLCDM model 
Q m (z) approaches 1 as soon as matter starts dominating, and this 
happens already for small values of z. From such z on, the value of 
/ becomes practically independent of 7 so that higher redshift data 
do not help constraining it anymore. If instead the model is more 
complex, e.g. flk is allowed to vary away from 0, then £l m (z) ~ 1 
at larger z so that increasing z max improves more the FoM. Some- 
thing similar also likely happens for fi, but it is more difficult to 
illustrate it since we only have numerical solutions for / in this 
case. Note also that for /i s there is more difference between curved 
and flat models than for 7. 

This and the fact that the FoM grows rapidly when increasing 
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Figure 13. Dependence of the FoM on the upper redshift limit of the survey z max (fixing the lower redshift limit z m j n to its minimum 0.7), computed for 
the growth parameters 7 (left panel) and fi s (right plot) and normalised to the FoM at Zmax = 1.0. The symbols have the same meaning as in Fig. ll2l but here 
the green solid line joins symbols corresponding to the most general model, curved CPL. 



Zmax encourages the effort to reach the highest possible maximum 
limiting redshift in future surveys. 

Finally, we test the dependence of the growth FoM on the red- 
shift interval covered by the survey. We plot in Figure [T4"1contours 
of constant FoM (in particular, FoM = 10%, 25%, 50% and 75% 
of the maximum FoM reached in each case) as functions of z m in 
and Zmax, where z min varies between 0.5 and 1.1 and z ma x be- 
tween 0.9 and 2.0. We concentrate here only on CPL dark energy, 
since this is the most complex and complete of our set of models. 
We first note that for the fj, s parameterisation there is stronger de- 
pendence on curvature of the FoM than for the 7 parameterisation 
(compare solid and dashed lines). We also note that the contours for 
7 and /is have different shapes. Let us consider a specific example. 
If we look at z m ;„ ~ 0.75, we need a z ma x > 1.6 to reach 50% 
of the maximum FoM, for both parameterisations. For other val- 
ues of Zmin, the two parameterisations give different indications. 
For Zmin < 0.75 (> 0.75), the minimum value of z max required 
in order to obtain 50% of the maximum achievable FoM is larger 
(smaller) for 7 than for fj, s . This study is warning us against opti- 
mising an experiment on the base of one parameterisation only. In 
the absence of a clear preference for one parameterisation over the 
other, an operative way to proceed in this particular case would be 
to choose the most conservative z m in and z max limits, for which 
the desired FoM is achieved in both parameterisations. 

To end this section on figures of merit, we list in Table [5] and 
[6] the FoMs, reached using the maximum available redshift range 



0.7 ^ z ^ 2.0 for all combinations of cosmological models, both 
for the 7 and the /i s parameterisation. We also look at how the FoM 
is improved when adding data from a low-z galaxy survey (as may 
be BOSS) and further add Planck data. In square brackets we show 
how the FoMs are degraded in the pessimi stic case of t he galaxy 
number density being half that forecasted in lGeachetalJj2008h . 



7 DISTINGUISHING GENERAL RELATIVITY FROM 
MODIFIED GRAVITY MODELS 

We finally turn to trying and answering the question "is the Euclid 
spectroscopic survey able to distinguish between GR and modified 
gravity?" To do this we use model selection tools (aimed precisely 
at telling how strongly a set of data prefers a model over other mod- 
el s) from Bayesian stat istics. In particular, we apply the method 
of iHeavens et alj J2007t) . where the Fisher matrix approach is gen- 
eralised to the context of model selection. The Bayesian tool for 
model selection is the Bayes' factor B, defined as the ratio of prob- 
abilities of model M' to model M, given the same data D, inde- 
pendently of the values assumed by the model parameters 6' or 6: 

_ P (M'\D) _ p(M')fd8'p(D\8',M')p(6'\M') 

p(M\D) p(M) J ddp(D\6,M)p(0\M) ( ' 

where p(M\D) is the probability of model M given the data D, 
p(M) is the prior probability of model M (i.e. the probability that 
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Figure 14. Dependence of the growth FoM on the redshift range for the two different growth parameterisations considered: contours of FoM = 10%, 25%, 
50% and 75% of the maximum FoM reached in each case (corresponding to the use of the maximum redshift range 0.7 ^ z ^ 2.0) are shown as a function of 
the lowest (z m i n ) and highest (z max ) redshift bins used. Here, a CPL model (dark energy with variable equation of state w = wq + w a (1 — a)) is assumed. 
Solid (dashed) lines correspond to a cosmology with = (Qk 0) Left panel: 7 parameterisation. Right panel: fi s parameterisation. 



Table 5. Figures of Merit for O m and 7. Both the case of fixing = (flat space) and allowing it to vary (curved space) were listed. All figures in square 
brackets represent the case of the galaxy number density being halved. The addition of other surveys at lower redshift was considered for all models, while the 
effect of adding Planck was computed only for one representative case, i.e. that of the most complex model (curved CPL). 



7 


Euclid 




+ low-z data 




+ Planck 




flat space 


curved space 


flat space 


curved space 


curved space 


qLCDM 


545 [361] 


209 [135] 


561 [376] 


221 [137] 




wCDM 


217 [146] 


74 [48] 


225 [153] 


75 [49] 




CPL 


35 [23] 


30 [20] 


35 [23] 


30 [20] 


141 [140] 



model M is true before the experiment is done), which is unknown, 
and we will assume p(M) = p(M'); p(D\9, M) is what is usually 
called the likelihood function (i.e. the probability that the data are 
true given the model M with parameters 9); p(9\M) is the prior 
probability of the model parameters 9 (i.e. the probability distribu- 
tion that one believes the model parameters have before the exper- 
iment is done). As in lHeavens et alj j2007f) . M' is here a dark en- 
ergy model well described by the CPL parameterisation, while M 
is a modified gravity model whose background expansion can still 
be described by the same H(z) of the CPL parameterisation, (not 
necessarily by the same fiducial parameters) but now the growth 
index 7 is fixed to a fiducial value different from the GR value of 
0.545. Our fiducial M' cosmol ogy is that of Eqs. {5). For our case 
the forecasted Bayes' factor is dHeavens et al.ll2007l) 



(B) = (2tv) 



-1/2 



VdetF 
Vdet F' 



exp 



A 7 (9) 



where A7 corresponds to a uniform prior on 7, F' and F are the 
Fisher matrices for GR and the modified gravity model, respec- 
tively, and 89 a represent the shifts of the parameters of M' due to 



a shift in 7 from the fiducial value 0.545: 

89 a = -(F'- 1 ) ap Gp{ 1 - lG R) (10) 

for a corresponding to all parameters but 7, while 89 a — 
(7 — 7gb ) for a associated to the 7 parameter. Here Gp is the 
vector drawn from F by extracting the column corresponding to 
the parameter 7 and the rows corresponding to all parameters ex- 
cept 7. From the above we see that (B) depends on the offset of 
7 with respect to the GR value, (7 — 'jgr), and on the prior on 7, 
A7. 

In Figure Q3] we show the dependence of the log of (B) 
on 1 7 — 7gr| for two different priors (represented by the solid 
and red dashed lines) in the case of our Euclid-like spectroscopic 
galaxy survey. The horizontal dotted lines describe values of In B 
which correspond to 'substantial' (1 < In B < 2.5), 'strong' 
(2.5 < lnB < 5) and 'decisive' (InB > 5) evidence in favour of 
one model with respect to t he other (bottom to top) according to Jef- 
freys' scale dJeffrevsll998l As can be seen from the figure, the Eu- 
clid spectroscopic survey alone will be able to substantially distin- 
guish between GR and a modified gravity model if 7— 7gr 0. 13, 
while it will be able to decisively distinguish between models if 
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Table 6. Figures of Merit for Q m and fj, s . Both the case of fixing = (flat space) and allowing it to vary (curved space) were listed. All figures in square 
brackets represent the case of the galaxy number density being halved. The addition of other surveys at lower redshift was considered for all models, while the 
effect of adding Planck was computed only for one representative case, i.e. that of the most complex model (curved CPL). 
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251 [165] 


94 [60] 
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Figure 15. Dependence of ln{B) on |7 — *)gr I f° r different priors in the Euclid spectroscopic galaxy survey. The solid black line corresponds to a uniform 
prior distribution for 7 with A7 = 1.0, the red dotted line to A7 = 0.7. The dotted red lines correspond to In B = 1, In B = 2.5 and In B = 5, delimiting 
the regions where evidence in favour of one model with respect to the other is 'substantial', 'strong' and 'decisive' according to Jeffreys' scale. The cusp 
corresponds to the case where B = 0, i.e. there is no evidence in favour of one model with respect to the other. This means that to the left of the cusp GR is 
favoured with respect to modified gravity models, while to its right modified gravity models are favoured. 



7 — 7GR 2? 0.2. We have also computed the evidence using the 
H parameterisation, with uniform prior distributions in the interval 
An a — 3 and A[i 3 — 5. For both priors it results that with Eu- 
clid spectroscopic data alone 'substantial' ('strong') evidence can 
be obtained in favour of a modified gravity if the latter has fi s > 1 
(/i s > 1.4). The addition of the weak lensing data from the Euclid 
photometric survey is expected to improve these results consider- 
ably. 



8 CONCLUSIONS 

In this paper we have investigated how strongly the Euclid galaxy 
spectroscopic survey in the current reference configuration can con- 
strain the growth of structure and consequently how well it can dif- 
ferentiate a GR cosmology from alternatives to it. 

We have found that we can reach precisions between 1.3% 
and 4.4% in the measurement of / erg depending on the redshift 
bin, where the highest precision is reached for z ~ 1.0. 

Comparing the Euclid spectroscopic survey with other ongo- 



ing and future galaxy redshift surveys we note that Euclid will reach 
the highest precision in the growth rate measurement. Euclid will 
be perfectly complementary to BOSS and BigBOSS: the three sur- 
veys together will allow to cover an extremely large redshift range: 
O.K z < 3.5. 

This precision in / as translates into a precision in the mea- 
surement of the growth index 7 which depends on the specific back- 
ground cosmology adopted. We have obtained marginalised errors 
on 7 — f2 m (or 7 — wo) between 5% and 10%. The parameterisation 
of the growth rate / we have adopted is / = ttZi (for curved space, 
/ = Q^i + (7 — 4/7)f2fc), where a departure from GR is repre- 
sented by a deviation of 7 from 0.545. We have considered nested 
background models: qLCDM (a model with constant w — —0.95), 
wCDM and CPL, both flat and curved. 

We have compared the relative gain in growth FoM (quantify- 
ing the precision in the joint measurement of Q. m and 7) for two 
different growth para meterisations, being the already mention ed 7 
and the parameter fj, a dPogosian et al. 2010: iSong etafll201ll) . We 
have found that when increasing the survey area (and correspond- 
ingly reducing the galaxy number density, having fixed the total 
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observing time) the FoM grows linearly. Moreover, this growth is 
quite mild. 

The relative gain in FoM when increasing z max is large, both 
using fi s and 7. We also have noted that (for curved models) the 
relative improvement increasing z max is approximately linear, en- 
couraging the effort to reach the maximum possible limiting red- 
shift. We then have examined the dependence of the FoM on the 
redshift interval covered by the survey. From Fig. |T4]it is possible 
to notice that to reach a desired FoM improvement one needs lower 
Zmin (or higher z max ) when considering /i s than when choosing 7. 
This warns us against relying on one parameterisation only when 
optimising an experiment. 

Finally we have forecasted the Bayesian evidence for Euclid 
and found that the spectroscopic survey alone will be able to sub- 
stantially (decisively) distinguish between GR and modified gravity 
models having 7 — 7gr ^ 0.13 (^ 0.2). This result is expected 
to improve even further when adding data from the photometric 
survey of Euclid, which improve not iceably the precisio n in the 
measurement of 7 (see e.g. Fig.2.5 of lLaureiis et alj|201 lb . 
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APPENDIX A: GROWTH OF STRUCTURES 

In a Friedmann-Lemaitre-Robertson- Walker (FLRW) universe, the 
growth of linear perturbations is described by the perturbed Ein- 
stein equations and conservation of the energy-momentum tensor 
equations. If we consider scalar perturbations about a FLRW back- 
ground in the Newtonian gauge, in Fourier space and in the limit of 
small scales, an equation for the growth rate / can be derived: 



df 
din a 



+ / 2 + 



2 + 



d\nH 
din a 



/ = ~n» 



(Al) 



A solution for / has now to be found and depends on the form of 
H. 



Al Parameterisations of growth in GR and deviations from 
it: a constant growth index 

A 1.1 Flat LCDM and quintessence 

If we assume a flat model with a simple cosmological constant or a 
dark energy with (quasi) constant equation of state w, with Fried- 



mann equation H = HQ[Q. m a + (1 — Q r , 



-3(1+™) 



], then it 



is possible to find a good solution of Eq. JAU by assuming 

f = Q m {aY', (A2) 

w here the growth index 7 is a constant. This was first proposed 
bv lPeebles! ( 1980) for a matter-dominated universe at z = 0, with 
7 ~ . 6. A b etter approximation for the same m odel was found 
bv lFrvl dl985b and lLightman & SchechteJ dl990|) . with 7 = 4/7. 
For the case of LCDM, 7 = 0.545 is a very good approxima- 
tion, at least in the redsh ift range of interest here. As shown in 
Wan g & Steinhardtl (l998), if we consider, instead of a standard 
cosmological constant, a quintessence field with slowly varying w, 
then 7 = 0.545 is still a good approximation, since the correction 
term has a weak dependence on w if this is close enough to —1. 

For larger redshift ranges, from Zqmb to today, and demand- 
ing a very high precision, a different solution for LCDM was 



proposed by Ishak & Dossett j2009h . For models wit h non-slowly 
varying w, iLinder & Cahnl d2007t) and iLindeil {2005 ) propose dif- 
ferent parameterisations. 

We restrict ourselves to models with slowly varying w and 
decide not to use the large redshift range parameterisation since we 
do not require an accuracy that high. Instead, we are interested in 
parameterisations for models where curvature is present, since this 
parameter has been shown to be degenerate with dark energy and 
we wish to explore it in a consistent way. 



A 1.2 Generalisation to curved LCDM and quintessence 

In the case of a curv ed LCDM mod el, t he first attempts at m odelling 
/ were proposed by Martel ( 1991) and La hav et al.ldl99ll) . 

More recently, in Gong et al. (2009]), two possible approaches 
are suggested. The first (less accurate) approach consists in still 
assuming f(z) — £l m (z) 7 , but changing the form of 7. The sec- 
ond consists in taking f(z) — fi 7 (z) 7 + otClk(z) (with Qk(z) = 
p k {z)/p{z)). 

In the second case, which is found to be more accurate, the 
solution is 

f(z) = iT m + (7 - 4/7)Q k , 7 = 0.545 . (A3) 

We use the parameterisation JA3) for our forecasts. 
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A1.3 Modified gravity: the case ofDGP, f(R) gravity and 
scalar-tensor theories 

A parameterisation as that of Eq. JA2) or jA3) is valid also for 
some modified gravity models (see Amendola & Tsuiikawa] |2010l 
for a complete review of viable alternatives to GR), when taking a 

different fiducial 7. 

Let us first take the DGP model. For flat space. lLinder & Cahnl 
J2007h showed that the growth rate / is still well par ameterised 
by Eq . jA2b with a constant 7. IWeil d2008h and iGong et al.l 
( 2009) derived the expe cted value of 7 b y solving the modi- 
fied g rowth equations of iLue et all J2004I) ; iKovama & Maartens 
3q): "1 = 11/16. A constant 7 still fits the curve well enough 
Linder & Cahnl 120071) . If we include curvature, we can proceed, 
following IGong et al. U2009h . as for curved quintessence, e.g. pa- 
rameterising f(z) as in Eq. jA3\ , and we find that the best fitting 7 
is again 7 = 11/16. 

Instead of adding extra dimensions, it is possible to modify 
ravity in four dimensions. This is done e . g. in f(R) theorie s 
Capozzielld 120021 : ICapozziello et alj 120031 : ICarroll et al. I l200i) . 
where the Einstein-Hilbert Lagrangian density R — 2A is mod- 
ified into a different function, R + f(R ). The function f(R) is 
strong ly constrained by local gravity tests dAmendola & Tsu iikawa 
2008), nevertheless the phenomenology of the growth of structure 
for these models is still very rich . Their growth index 7 is in gen- 
eral t i me- and scale-dependent jTsuii kawa 2007; Gannou il et al.l 
20091: iTsuiikawa et alj |2009|; iMotohashi et alj|20ld) . As found in 
Tsuiikawa et alj d2009f ) and lGannouii et al.l 1 20091) . the importance 



of scale-dependent dispersion of 7 depends on the particular f(R) 
model chosen and on its parameters. For some particular values 
of the parameters, there can be scale-independent growth. We do 
not consider scale dependence in this work, although it would be 
very interesting to ex amine this feature, as it might be a smoking 
gun for /(i?) models (Pogo sian et alj |2010). Time dependence can 
arise in general. Forecasts on 7(a) for Euclid - like g alaxy redshift 
surveys have been computed in lDi Porto et alj d201 ll) and we refer 
the reader to this paper for a full analysis of the topic. In this work 
we concentrate on a constant 7, which somehow corresponds to a 
7(a) averaged over the redshift range of the survey, weighted by 
the error on / in each redshift bin. Given that viable f(R) models 
show a lower growth index 7(2 = 0) ~ 0.4 with respect to GR, be- 
coming even smaller at higher z, if one was to detect an unusually 
low averaged constant 7, this would point to f(R) models. The use 
of a time and possibly scale dependent 7 could then in principle 

allow to distinguish among different f(R) models. 

S calar-tensor theories dAmendola! 1 19991: Uzar]|l999t IChibal 



1 19991 : iBartolo & Pietronil 120001; IPerrotta et al.l I2OO0J) are a gen 
eralisation of f(R) models, where the Lagrangian density is 
1/2 R) - 1/2 C(0)(V0) • Given the generality of these the- 
ories, the phenomenology of their growth of structure is very rich. 
The growth of matt er perturbations in some of t h ese models has 
been s t udied e.g. inlDi Porto & Amend ola (2008); Gannou ii et al.l 
(2009); iKobavashil fcOld) . In general 7 depends on time (and 
scale), but t here are also spe cific cases where 7 ~ const for red- 
shifts < 2 dKobavashill20Ta) . 



stead of some already approximated version of them (e.g for small 
scales). This approach may turn useful when combining constraints 
from different observational tools, which might need different ap- 
proximations. This way, the parameters would always have a physi- 
cal meaning and one would a void difficulties in inte rpreting results, 
as pointed out very neatly in lPogosian et alj d2010T>. It i s prec isely 
the parameterisation (fi, 77) proposed in Pogosia n et al.l J2OI0I) that 
we use and compare to 7. The Einstein equations defining the pa- 
rameterisation are: 

$ 



-AtyGno, fi(a,k)p5, 



r)(a, fc) , (A4) 



where and $ are the metric perturbations in the Newtonian 
gauge: 



ds 



-a [(1 + 2*)dr 2 



2§)dx 



(A5) 



Since the galaxy power spectrum is sensitive to \l/ and not to the 
anisotropic stress (which depends on the difference between the 
two potentials $ and <&), the parameter /j, is sufficient for our work 
and we do not use 77 at all. 

For simplicity, we de cide to assume again s cale-independence: 
(i — fi(a) (although see Pogo sian et alj|201fl and Se c. IA1.3I for 
reason s to keep scale-dependence). We model fj, as in ISong et alj 
d201ll): 



H(a) = 1 + fi s a , 



(A6) 



which is motivated by DGP and reduces to GR for /j, s = 0. So, 
using again the function / defined previously, we obtain 

jlL + f + (2 + ^)f-ln m (l + Vsa) = 0. (A7) 
a In a \ a In a J 2 

We can solve this equation numerically, by imposing the ini- 
tial condition at an initial redshift of matter domination z m d'- 
f (zmd) = QlzrndY 1 (w here 7 can be determined from Eq. (32) 
of IPogosian etail d20ldV ). 



APPENDIX B: PREDICTIONS USING THE /j 
PARAMETERISATION 

Here we present forecasts on the marginalised errors on cosmo- 
logical parameters when, instead of using 7, the alternative fj, pa- 
rameterisation is used (see Appendix IA2t . As we can see from 
Figs. IB 1 1 and IB2I with Euclid spectroscopic data alone the abso- 
lute marginalised error on /i s and Q m or wo will be ~ 2. Adding 
Planck will improve the measurement of /i s by more than 50%, 
reducing the error to ~ 0.7 (see Figs. |B3| and lB4l ). 

This paper has been typeset from a TgX/ LTj-X file prepared by the 
author. 



A2 Parameterisations of growth in GR and deviations from 
it: physical parameters 

Instead of parameterising the deviati on from GR with 7, w e can 
use a different approach, proposed bv lAmendola et al1 (2008). This 
consists in directly parameterising the full Einstein equations in- 
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Figure Bl. Forecasted marginalised errors on f2 m and fj, s using the Eu- 
clid spectroscopic survey. Solid (black) lines correspond to qLCDM, dashed 
(red) lines to wCDM and dot-dashed (cyan) lines to CPL. Thick lines rep- 
resent models with Qj. = 0, while thin lines indicate curved models. 
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Figure B4. Comparison of marginalised errors on uiq and /i s in curved 
CPL, obtained using Euclid only and joint Euclid and Planck data. Dot- 
dashed (cyan) contours correspond to error ellipses using Euclid data only, 
while solid (purple) contours show joint constraints from Euclid and Planck. 
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Figure B2. Forecasted marginalised errors on turj and fj, s using the Euclid 
spectroscopic survey. Dashed (red) lines correspond to wCDM and dot- 
dashed (cyan) lines to CPL. Thick lines represent models with Qj. = 0, 
while thin lines represent curved models. 
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Figure B3. Comparison of marginalised errors on Q m and fi s in curved 
CPL, obtained using Euclid only and joint Euclid and Planck data. Dot- 
dashed (cyan) contours correspond to error ellipses using Euclid data only, 
while solid (purple) contours show joint constraints from Euclid and Planck. 
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